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Previously published data on inclusive neutral and charged current e^p and e~p deep inelastic scattering from 
HERA have been combined. The new, model-independent method of combining the measured cross sections 
takes full account of correlated systematics in a coherent way, leading to significantly reduced uncertainties in 
the combined measurement across the {x,Q^) plane. The combined cross section data have been used as the sole 
input for a new next-to-leading order QCD fit to extract the parton distribution functions (PDFs) of the proton. 
The PDFs extracted have greatly reduced experimental uncertainties, compared to separate QCD analyses on 
data from HI and ZEUS. Model uncertainties, including those arising from the parameterisation dependence, 
have also been carefully considered. The resulting HERA PDFs have impressive precision. 



1. Introduction 

The deep inelastic lepton-nucleon scattering 
(DIS) process, in which a point-Uke lepton probes 
a single, initial state hadron, is an ideal environ- 
ment in which to study quantum chromodynam- 
ics (QCD), and constrain the parton distribution 
functions (PDFs) of the proton. 

The HERA collider ceased running in June 
2007, following 15 years of successful operation. 
During HERA I running (1992-2000), the HI and 
ZEUS experiments each collected approximately 
100 pb^^ of e+p and ^ 15 pb^^ of e~p data. This 
has allowed precision measurements of inclusive 
neutral (NC) and charged (CC) current DIS cross 
sections, which have already proved fundamen- 
tal to the rapid development in understanding of 
QCD and the structure of the proton. 

In these proceedings, a joint HI and ZEUS 
analysis is presented, in which previously pub- 
lished NC and CC inclusive DIS cross section 
measurements from the two collider experiments 
are combined. The combination method [l] 
uses an iterative minimisation, which care- 
fully takes into account the correlations within 
the data that result from different sources of un- 
certainty. The key assumption is that the mea- 
surements from HI and ZEUS represent a com- 
mon truth. Thus, forcing them to agree results 
in a strong constraint which cross-calibrates the 



measurements, resulting in significantly reduced 
overall uncertainties. The combined HERA data 
have subsequently been used as the sole input to a 
next-to-leading (NLO) QCD fit to determine the 
proton PDFs. This analysis is also discussed in 
the present contribution. 

The results on the combined HERA data, as 
presented here, were first released to the LP Con- 
ference in 2007 [2], and the NLO QCD analysis 
was first presented at DIS 2008 3 . Both anal- 
yses should be seen as part of a mid-term strat- 
egy. They will be followed by future data combi- 
nations, including even more accurate data from 
both HERA I and HERA II, and by further QCD 
analyses to extract the proton PDFs. 

2. HERA physics and kinematics 

Lepton-proton DIS can proceed either via the 
neutral current (NC) interaction (through the ex- 
change of a 7* or Z°), or via the charged cur- 
rent (CC) interaction (through the exchange of a 
W^). The kinematics of lepton-proton DIS are 
described in terms of the Bjorken scaling vari- 
able, X, the negative invariant mass squared of 
the exchanged vector boson, Q^, and the fraction 
of energy transferred from the lepton to the pro- 
ton system (in the rest frame of the proton), y. 
The three quantities are related by = s ■ x ■ y, 
where s is the centre-of-mass energy squared. 
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At leading order (LO) in the electroweak inter- 
action, the reduced cross section, ar{x,Q^), for 
the e^p NC DIS process can be expressed in terms 
of proton structure functions, 

,2 
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1 dV^ 



NC 



27ra2 Y+ dxdQ'^ ' 

where a is the fine structure constant and Y± — 
1 ± (1 — y)^ describes the helicity dependence of 
the electroweak interaction. The structure func- 
tions, which depend on (x, Q^), are directly re- 
lated to PDFs of the proton, and their de- 
pendence is predicted by perturbative QCD. In 
particular, F2 and xF^ depend directly on the 
quark distributions. For < iqOO GeV^, F2 
dominates the ep scattering cross section and for 
X < 10"^, F2 itself is dominated by sea quarks 
while the dependence is driven by gluon ra- 
diation. F2 has been measured over 4 orders of 
magnitude in {x, Q^), to a precision of ~ 2 — 3% 
at HERA I. Therefore, HERA data have already 
provided vital information on the sea-quarks and 
gluon at low x. At high > Al'^, the contri- 
bution from xF^ becomes increasingly significant 
and gives information on the valence quarks. 
At LO, the CC cross sections are given by, 



, + c+{l^vf{d + s)\ 
^ + c+{l-yf{d+-s)\ 



dxd02 2^(g2 + M^Y 

dxdQ^ 2^(g2 + M^Y 

SO that a measurement of the e^p and e^p cross 
sections provides information on the d- and u- 
valence quarks, respectively, thereby allowing the 
separation of flavour. 

3. Combination of HERA DIS data 

In this section, the joint HI and ZEUS com- 
bined data analysis is described. The goal of the 
study is to obtain DIS cross sections of best possi- 
ble accuracy in order produce precise extractions 
of the proton PDFs. A new NLO QCD analysis, 
in which the results of the combination are used 
as the sole input, is described in Sec. ID 



3.1. Combination Method 

The averaging procedure uses the Lagrange 
Multiplier (or Hessian) method [4 . The only the- 
oretical input to the combination is that there is 
a true value of the cross section for each process, 
at each (a;,^^) value |2 5j. The correlated sys- 
tematic uncertainties are floated coherently such 
that each experiment calibrates the other one. 
This results in a significant reduction of the cor- 
related systematic uncertainty over much of the 
kinematic plane. 

In the combination procedure, the following 
probability distribution of a measurement quan- 
tity, Af, represented as a function, is min- 
imised: 



Xea; 
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Here, A'P are the measured central values, and 
Si are the statistical and uncorrelated systematic 
uncertainties of the quantity AI. The are 
their true values; aj are parameters for the j 
sources of systematic uncertainty and dAP/daj 
denotes the sensitivity of point i to source j. For 
the cross section measurements, the index i la- 
bels a particular measurement at a given {x, Q^). 
Equation [T] represents the correlated probability 
distribution function for the quantity M*'*''"'' and 
for the systematic uncertainties aj. 

The defined in Eq. [T] has, by construc- 
tion, a minimum = for Af*'*'"" = Af* and 
aj ~ 0. The total uncertainty for M*'*"^^" de- 
termined from the formal minimisation of Eq. [T] 
is equal to the quadratic sum of the statistical 
and systematic uncertainties. The covariance ma- 
trix cov(Af''*''"°, Af'^'*'"'') quantifies the correla- 
tion between experimental points. 

In the analysis of more than one data set, a 
total function, Xtot^ constructed from the 
sum of the functions for each data set. The 
Xtot function can be minimised with respect to 
j^^i,truc g^j^j Q,^.. ^j-^jg minimisation corresponds to 
a generalisation of the averaging procedure which 
takes account of correlations between different 
data sets. 
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data set 






range (GeV^) 


L. (pD ) 


comment 


ref. 


HI NC min. bias 


97 


0.00008 


0.02 


3.5 


12 


1.8 






301 GeV 


m 


HI NC low Q2 


96-97 


0.000161 


0.20 


12 


150 


17.9 




Vi = 


301 GeV 


m 


HI NC 


94-97 


0.0032 


0.65 


150 


30000 


35.6 


e+p 




301 GeV 


m 


HI CC 


94-97 


0.013 


0.40 


300 


15000 


35.6 


e+p 


V~s = 


301 GeV 


m 


HI NC 


98-99 


0.0032 


0.65 


150 


30000 


16.4 


e^p 




319 GeV 


m 


HI CC 


98-99 


0.013 


0.40 


300 


15000 


16.4 


e~p 




319 GeV 


m 


HI NC 


99-00 


0.00131 


0.65 


100 


30000 


65.2 


e+p 




319 GeV 


m 


HI CC 


99-00 


0.013 


0.40 


300 


15000 


65.2 


e'^p 




319 GeV 


m 


ZEUS NC 


96-97 


0.00006 


0.65 


2.7 


30000 


30.0 


e^p 




301 GeV 


m 


ZEUS CC 


94-97 


0.015 


0.42 


280 


17000 


47.7 


e+p 




301 GeV 


m 


ZEUS NC 


98-99 


0.005 


0.65 


200 


30000 


15.9 


e^p 




319 GeV 


m 


ZEUS CC 


98-99 


0.015 


0.42 


280 


30000 


16.4 


e^p 




319 GeV 


m 


ZEUS NC 


99-00 


0.005 


0.65 


200 


30000 


63.2 


e+p 




319 GeV 


m 


ZEUS CC 


99-00 


0.008 


0.42 


280 


17000 


60.9 


e+p 




319 GeV 


m 



Table 1 

The HI and ZEUS datasets used in the combination procedure. The integrated luminosity [C) and the 
kinematic range in {x, Q^) are given. Note that a re-analysis of the HI luminosity measurement of the 
special minimum bias run in 1997 has lead to a change in integrated luminosity. This has resulted in an 
upward shift of the HI minimum bias data by 3.4%, which is taken into account in the combination. 



The function of Eq.[l]is most suitable for abso- 
lute or additive uncertainties, i.e. those which do 
not depend on the central value of the measure- 
ment. However, for cross sections, many uncer- 
tainties are proportional to the central value (so- 
called multiplicative uncertainties). This propor- 
tionality can be approximated by a linear depen- 
dence. In such cases, the data combination using 
Eq. [1] will introduce a bias towards lower cross 
sections since the measurements with smaller cen- 
tral values will have smaller absolute uncertain- 
ties. An improved can be defined by re- 
placing, in Eq. ^ 5i ^^jp^Si and f^^j ^ 

— '^J which translates the multiplicative 
uncertainties for each measurement to the ab- 
solute ones. Unlike Eq. [l] however, this is 
not a simple quadratic function with respect to 
{M*'*''"'^},{aj}. Therefore, the minimum is found 
by an iterative procedure: first, Eq. [1] is used 
to find an initial approximation for {JVP'*^^'^}, 
then the uncertainties are re-calculated using 
di Mi — finally the minimisation is 
repeated. For the HERA data averaging, conver- 
gence was observed after two iterations. 



3.2. Data input and treatment 

The data used for the combination consist of 
the published double differential NC and CC cross 
sections from HI and ZEUS, taken in the years 
1994-2000, and are listed in Tab. [2 During this 
period, HERA operated with an electron beam 
energy, E^, of 27.5 GeV and a proton beam en- 
ergy, Ep, of 820 GeV (until 1997) and 920 GeV 
(from 1998 onwards). The measurements span 
the kinematic region 1.5 < < 30000 GeV^ 
and 6.5 x 10~^ < a; < 0.65, and are the most 
precise data published by the HI and ZEUS col- 
laborations to date0. 

The double differential cross section measure- 
ments are published with their statistical and 
systematic uncertainties. The statistical uncer- 
tainties are uncorrelated between different data 
points. The systematics are classified into three 
sub-groups: (i) point-to-point uncorrelated sys- 
tematics, (ii) point-to-point correlated (e.g. en- 



^ Note that there are also data available for < 1 GeV'^, 
both from shifted vertex operation and from ZEUS using 
a dedicated detector near the beam pipe. These data have 
not been considered here, but will be included in subse- 
quent combined data analyses. 
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HERA I e*p Neutral Current Scattering - HI and ZEUS 




Q-/GeV- 



Figure 1. Deep inelastic neutral current e^p scat- 
tering cross section data from the HERA I data- 
taking period as obtained by combining the pub- 
lished HI and ZEUS measurements. The curves 
are NLO QCD fits as performed by HI and ZEUS 
to their own data. 
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Figure 2. Neutral current e~^p cross section mea- 
surements for three selected x bins as a function of 
Q^. The HI (open points), ZEUS (open squares) 
are compared to the HI and ZEUS combined data 
(solid points). Measurements from the individual 
experiments have been shifted for clarity. The 
error bars indicate the total uncertainty. The 
curves show the predictions of NLO QCD fits per- 
formed by HI and ZEUS, on their own data. 



ergy scale calibrations), (iii) an overall normali- 
sation uncertainty of various data sets. Sources 
of point-to-point correlated uncertainties are of- 
ten common for NC and CC cross section mea- 
surements and sometimes can be considered to 
be correlated for different data sets of the same 
experiment. They are treated as independent be- 
tween HI and ZEUS, since uncertainties of beam 
energies are negligible. Similarly, the normalisa- 
tion uncertainties are correlated for all cross sec- 
tion measurements by a given experiment from a 
common data-taking period. 

All the NC and CC cross section data from 
HI and ZEUS are combined in one simultaneous 
minimisation. Therefore, resulting shifts of corre- 
lated systematic uncertainties and global normal- 
isations propagate coherently to both NC and CC 
data. 



3.2.1. Extrapolation to common x — 

Prior to the combination, the HI and ZEUS 
data were transformed to a common grid of 
{x,Q^) points using a simple interpolation: 

ep(X grid, W grid ) ~ rrth ( ^ n'^\ <^ep\X,y^ )■ 
"ep y^i^ ) 

The HI PDF parameterisation [16] of the double 
differential NC and CC cross sections was used to 
calculate the theoretical ratios. The sensitivity of 
the data combination to the choice of parameter- 
isation was checked using the ZEUS PDFs [T7] . 
The resulting correction factors were found to 
agree to within a few permille for the NC data 
and to within 2% for the CC data (i.e. in both 
cases, much less than the experimental uncertain- 
ties). 
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3.2.2. Centre-of-mass energy correction 

The data sets listed in Tab. [U include sam- 
ples from both Ep = 820 GeV and Ep = 920 GeV 
running. Since the NC and CC DIS cross sections 
depend weakly on the energy, via terms involving 
the inelasticity y, a choice must be as to whether 
to keep the two beam-energy samples separate, 
or to correct to a single, common proton beam 
energy. For the present combination, the latter 
choice was made such that the NC and CC data 
at Ep = 820 GeV were transformed to Ep ~ 920 
GeV and then combined with the genuine mea- 
surements at the higher proton beam energy. 

For the CC data , the correction was calculated 



as: 



'''(7(7,920 (•^J 



th^ep 
^CC,920 



•^C(7,82o('''' Q^) 



ep 
^CC,820 



For the NC data, the correction was performed 
additively: 



ep 

^NC,920 



+ Acr^pg2o(a;,Q^y92o, 2/820 ), 

where the correction term is the difference be- 
tween HI PDF [in] based predictions of the DIS 
cross sections, with 1/920 — / {4:xEe920) and 
2/820 = /{4:xEe820). The corrections were 
found to only be sizable at large y. To estimate 
the uncertainty on the combined data due to this 
procedure, another average was calculated assum- 
ing = as an extreme assumption. The dif- 
ference between these results and those using the 
standard choice were at the permille level across 
most of the kinematic plane, and reached 5% 
for only a very few points at y > 0.6. This uncer- 
tainty was added in quadrature to the combined 
data. 

3.2.3. Additive vs multiplicative errors 

As mentioned in Sec. 13. 11 uncertainties on cross 
section measurements may be additive or multi- 
plicative. Within the community it is generally 
agreed that normalisation uncertainties are mul- 
tiplicative. However, for the other systematic un- 
certainties the situation is less clear. To assess the 
sensitivity of the average HERA data set to this 
issue, various different treatments of the system- 
atic uncertainties were considered. The extreme 



assumptions treat all uncertainties as multiplica- 
tive, or all as additive, apart from the normal- 
isation uncertainties. Therefore, an additional 
systematic uncertainty was estimated, based on 
the difference between these two error treatments. 
The typical size of this uncertainty was < 1% for 
the low data, reaching 1 — 1.5% at larger Q^. 

3.2.4. Correlations between experiments 

The HI and ZEUS collaborations use similar 
methods to reconstruct the event kinematics, em- 
ploy similar techniques for the detector calibra- 
tion, use common Monte Carlo simulation mod- 
els for the hadronic final state simulation as well 
as for photoproduction background subtraction. 
This similarity of approaches and techniques may 
lead to correlations between HI and ZEUS mea- 
surements. 

A detailed investigation has shown that the re- 
sults of the combination are rather insensitive to 
the assumptions on correlations between the two 
experiments. The largest effect on the average de- 
rives from differing assumptions on the photopro- 
duction background (a 1 — 2% change at y > 0.6 
for low < 20 GeV'^) and on the hadronic en- 
ergy calibration (1% at low y < 0.02). For these 
sources the measurements rely more on the sim- 
ulation of the hadronic final state which is simi- 
lar for the two experiments. These variations are 
therefore introduced as additional point-to-point 
correlated systematic sources of uncertainty on 
the averaged cross sections. 

3.3. Results 

In the minimisation procedure, 1153 individ- 
ual NC and CC measurements were averaged to 
584 unique points. This yielded a good qual- 
ity of fit with the x^/dof = 510/599. The dis- 
tribution of pulls did not show any significant 
tension across the kinematic plane. A total of 
43 sources of correlated systematic uncertainty 
were considered in this analysis. In the com- 
bined data, almost all systematic uncertainties 
were reduced, with the most significant reduc- 
tion observed for the HI backward calorimeter 
energy scale (by a factor of 3) , and the ZEUS un- 
certainty in modelling the forward hadronic en- 
ergy flow (by a factor of 4). The errors which re- 
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suit from the combination process (centre-of-mass 
energy correction, additive versus multiplicative 
error treatment, correlations in the background 
subtraction and on hadronic energy calibration) 
were introduced as additional point-to-point cor- 
related systematic sources of uncertainty on the 
average cross section. 

Figure [1] shows the complete set of combined 
e'^p NC data, spanning the entire x and range. 
The scaling violations are clearly visible at both 
low and high x. The data are compared to two 
previously pubhshed QCD fits performed by the 
HI and ZEUS collaborations on their own data. 
As expected, the fits provide an excellent descrip- 
tion of the HERA combined data. 

Figure [2] shows a close up, with a linear y scale, 
of three selected x bins, as a function of Q^. The 
combined data are compared to the individual 
measurements from HI and ZEUS (shifted for 
clarity). At low Q^, where the data are limited 
by systematic uncertainties, the improvement in 
the total uncertainty is visible. At higher Q^, 
the combined data have significantly reduced un- 
certainties, and exhibit far smaller fluctuations, 
which is driven by the increase in statistical ac- 
curacy, which dominates the measurement. 

The combined NC and CC HERA data have 
subsequently been included in a NLO QCD fit to 
extract the proton PDFs. Details of this analysis 
are described in the next section. 

4. QCD analysis of the combined data 

Previously, the HI and ZEUS collaborations 
have both used their own data in NLO QCD 
fits [16|17j . These data sets have very small statis- 
tical uncertainties, so that the contribution from 
systematics becomes dominant and consideration 
of point-to-point correlated systematic uncertain- 
ties is essential. The ZEUS analysis takes ac- 
count of correlated experimental systematics by 
the Offset method (see e.g. [H]), while HI uses 
the Hessian method [?] . In an attempt to improve 
the determination of the PDFs from HERA data, 
the combined HI and ZEUS measurements have 
been used as the sole input for a new, NLO QCD 
analysis. In Sec. 14. 1[ the QCD analysis, model 
assumptions and treatment of uncertainties are 



discussed. In Sec. 14.21 the results are presented. 

4.1. NLO QCD analysis 

For the QCD analysis presented here the pre- 
dictions for the structure functions were obtained 
by solving the DGLAP evolution equations at 
NLO, in the MS scheme, with the renormalisa- 
tion and factorisation scales taken to be . The 
DGLAP equations yield the PDFs at all values 
of provided they are input as functions of x 
at some input scale Qq. For this analysis, the in- 
put scale was chosen to be = 4 GeV^ The 
resulting PDFs were then convoluted with NLO 
coefficient functions to give the structure func- 
tions which enter into the expression for the cross 
sections. The choice of the heavy quark masses 
were TOc = 1-4 GeV and mt = 4.75 GeV. For this 
preliminary analysis, the heavy quark coefficient 
functions have been calculated in the zero-mass- 
variable-flavour-number scheme. The strong cou- 
pling constant was fixed at as{Mz) = 0.1176 [T8] . 

The fit was performed at leading twist. Since 
the HERA data have a minimum invariant mass 
squared of the hadronic system, W^, of VF^j„ = 
300 GeV^ and a maximum x of Xmax = 0.65, they 
are in a kinematic region where there is no sen- 
sitivity to target mass and large- a; higher twist 
contributions. However, a minimum cut of 
Qmin = 3.5 GeV^ has been imposed on the data 
included in the fit, in order to remain in the kine- 
matic regime where perturbative QCD should be 
applicable. 

4.1.1. Choice of parameter isat ion 

The PDFs were parameterised using the form: 

xf{x) = Ax^{l - x f{l + Dx + Ex^ + Fx^), 

and the number of parameters was chosen such 
that D, E and F were only varied if this resulted 
in a significant improvement to the x^. Other- 
wise, these parameters were set to zero. 

For the present analysi^, the following PDFs 

^Note that the choice of parameterisation described in 
Sec. 14.1.11 has been inspired by both the HI- and the 
ZEUS-style parameterisations. It attempts to combine 
the best features of both, in that it has fewer assump- 
tions concerning x{d — u) than the ZEUS-style and less 
model dependence than the Hl-style since it does not as- 
sume equality of all B-parameters. 
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model variation 


central value 


lower 


upper 


rric 


1.4 


1.35 


1.5 




4.75 


4.3 


5.0 


^ mm 


3.5 


2.5 


5.0 


Qo 


4.0 


2.0 


6.0 


fs 


0.33 


0.25 


0.40 


fc 


0.15 


0.12 


0.18 



Table 2 

Central values of input parameters and cuts, and 
the variations considered to evaluate model un- 
certainties. 



were parameterised: xuy, xd^, xg, xU — x{u + c) 
and xD = x{d + s) such that: 

= Au^x^""- {1 - x)^''^{l + Du^x + Eu^x'^) 

xd^x) = Ad^x^^- {1 - xf-'- 

xU{x) = %a;'^-(l - x)^- 

xD{x) = A^x^~{l - 

xg{x) = AgX^^{l - x)'^^ . 

The normalisation parameters, and Afi^ were 
constrained to satisfy the number sum rules and 
Ag was constrained to satisfy the momentum sum 
rule. The B parameters, Bu^ and Bd^ were set 
equal, i?„„ — Bd^ , as were the Bjj and B-j^, 
Bjj — Bjj, such that a single i?-parameter was 
used for the valence quarks and another, differ- 
ent single B-paramcter was used for the sea dis- 
tributions. A further constraint was provided by 
assuming that the strange and charm quark dis- 
tributions can be expressed as x dependent frac- 
tions, fs — 0.33 and fc = 0.15, of the d and u 
type sea, such that A{U) ^ A(D){1~ fs)/{l~ fc)- 
The value of fs = 0.33 was chosen to be consistent 
with determinations of this fraction from neutrino 
induced dimuon production. The charm fraction 
was set to be consistent with dynamic generation 
of charm from the start point of Qq — rri^ in 
a zero-mass-variable-flavour-number scheme. In 
total, there were 11 free parameters in the fit. 

4.1.2. Treatment of experimental data 

The HERA combined e'^p NC and CC inclusive 
data, described in Sec. [H have been used as the 



experimental input to the present QCD analysifl 
The consistency of the combined HERA data set, 
as well as its small systematic uncertainties, allow 
the experimental uncertainties on the PDFs to 
be calculated using a tolerance of Ax^ = 1- 
This is compared to previous global fits where 
increased tolerances of Ax^ = 50— 100 have been 
used in order to account for data inconsistencies. 

For the present QCD fit, the role of corre- 
lated systematics is no longer crucial since the 
uncertainties are relatively small. In particu- 
lar, this means that similar results are obtained 
whether the correlated systematic uncertainties 
are treated using the Offset or Hessian method in 
the QCD fit, or by simply combining statistical 
and systematic uncertainties in quadrature. For 
the central fit the choice was made to combine the 
43 systematic uncertainties (resulting from the 
separate HI and ZEUS data sets) in quadrature 
and to Offset the 4 sources of uncertainty which 
arise from the combination procedure. This was 
found to result in the most conservative uncer- 
tainty estimate on the extracted PDFs. 

4.1.3. Model uncertainties 

Despite the conservative procedure adopted, as 
mentioned above, the experimental uncertainties 
on the resulting PDFs were still found to be im- 
pressively small. Therefore, a thorough consid- 
eration of further uncertainties due to model as- 
sumptions was necessary. For the present anal- 
ysis, six sources of model uncertainty were con- 
sidered, as listed in Tab. [21 The parameters were 
varied up and down from the central value, and 
the differences in the resulting PDFs from the 
central fit were added in quadrature to the total 
experimental PDF uncertainty. 

Further cross checks were also performed. In 
particular, the dependence on the choice of pa- 
rameterisation was investigated by repeating the 
fit using the Hl-style [H] and ZEUS-style 
parameterisations. All fits gave acceptable 
and were found to be consistent with each other. 
Note, however, that the parameterisation chosen 

^Note that the version of the combined data used for 
the QCD analysis is slightly different to that described in 
Sec. |3] in that only data with y < 0.35 were transformed 
in E„. 
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Figure 3. HERA combined NC (left) and CC (right) reduced cross section data. The predictions of the 
HERAPDFO.l fit arc superimposed. The uncertainty bands iUustrated derive from both experimental 
and model sources. 




Figure 4. The HERA PDFs xuv, xdv, xS (total sea-quark distribution) and xg at (left) ~ 10 GeV^ and 
(right) = 10000 GeV^. Fractional uncertainty bands are shown beneath each PDF. The experimental 
and model uncertainties are shown separately. 
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Figure 5. Left: The HI and ZEUS PDFs obtained ft-om separate QCD analyses, each to their own data. 
Right: The HERAPDFO.l PDFs obtained from the analysis of the combined data set, as described here. 




Figure 6. The HERAPDFO.l PDFs compared to the CTEQ6.1M and MRSTOl fits at = 10 GeV^ 
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for the central fit (described in Sec. 14. ip gave the 
best weW as the most conservative experi- 

mental uncertainties. 

4.2. Results 

The NLO QCD analysis described above has 
been named the HERAPDFO.l fit. Figure [3] 
shows the results of the HERA fit superimposed 
on the combined high-Q^ NC and CC e^p data 
sets, showing the excellent description of the data 
by the HERA PDFs. 

Figure HI shows the HERAPDFO.l PDFs xu^, 
xdv, xS (total sea) and xg as a function of x at 
0^ = 10 GeV^ and = 10000 GeV^. Frac- 
tional uncertainty bands are shown beneath each 
plot, with the experimental and model uncertain- 
ties being shown separately. The variation of the 
strange fraction, fs, dominates the model uncer- 
tainty on the sea, while variations of Qq and Q^i„ 
dominate those on the gluon and valence quarks. 
Comparison of the results for = 10 GeV^ and 
0^ = 10000 GeV^ shows that, as increases, 
the PDF uncertainties become impressively small. 

The summary plots shown in Fig. [5] illustrate 
that the total uncertainty of the PDFs obtained 
from the HERA combined data set is much re- 
duced compared to the those of the PDFs ex- 
tracted from separate analyses of the HI and 
ZEUS data sets. The dramatic improvement is 
a result of the data combination. Figure [6] com- 
pares the HERAPDF 0.1 PDFs to those of the 
CTEQ6.1 [IS] andMRSTOl [ID] global fitfl The 
results indicate that the precision of the HER- 
APDFO.l PDFs for the low x sea and gluon is 
impressive. 

5. Summary 

A new, model-independent method of combin- 
ing cross section measurements has been pre- 
sented, in which a coherent treatment of system- 
atics results in a substantial reduction in the over- 
all uncertainties. The method has been demon- 
strated on the complete set of published HERA I 
NC and CC inclusive DIS data. 

^Note that the HERAPDFO.l uncertainty band represents 
a 68% confidence level, while the global fits show a 90% 
confidence level band. 



The combined HERA data have subsequently 
been included in a new NLO QCD analysis. The 
consistent treatment of systematic uncertainties 
in the combined data set ensures that the exper- 
imental uncertainties on the PDFs can be cal- 
culated without need for an increased toler- 
ance. This results in PDFs with greatly reduced 
experimental uncertainties compared to separate 
analyses of data from the HI and ZEUS ex- 
periments. Model uncertainties have also been 
carefully considered. The resulting PDFs (called 
HERAPDFO.l) have impressive precision. They 
are now publicly available in LHAPDF f2T (v5.6 
onwards). 

High statistics data from the HERA II running 
period are being analysed by the two collabora- 
tions and will be included in subsequent analyses 
devoted to precision determinations of the proton 
parton densities. 
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